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Abstract 

The marginal distribution of squeezed and rotated quadrature for two types of nonclassi- 
cal states of trapped ion - for squeezed and correlated states and for squeezed even and odd 
coherent states (squeezed Schrodinger cat states) is studied. The obtained marginal distri- 
Q\ ■ bution for the two types of states is shown to satisfy classical dynamical equation equivalent 

to standard quantum evolution equation for density matrix (wave function) derived in sym- 
plectic tomography scheme. 
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1 Introduction 



Recently, in |T[ it was shown that the steady state of trapped ion irradiated by bichromatic 
laser field is superposition of two coherent states [Q] which is even/odd coherent state introduced 
in H and interpreted in Jlj for large amplitudes of the superposition partners as Schrodinger cat 
states [|J. The theory of ion in the Paul trap was developed in [|, [7|, where the trapped ion 
was described by the model of quantum oscillator with periodically varying frequency. Gaussian 
packets, discrete modes and propagator of such oscillator have been obtained in Linear in 
position and momentum integrals of motion have been found for the oscillator in ||, |7|, [|. In fllCH , 
on the basis of analysis [fTI], the procedure was formulated to obtain the Wigner function of 
quantum system in terms of marginal distribution of rotated quadrature which may be measured 
by balanced homodyne detector. The scheme which was called optical tomography has been used 



experimentally In the symplectic tomography procedure has been suggested in which 
measuring the quantum states was proposed by means of measuring marginal distribution for 
squeezed, rotated and shifted quadrature. 
^ In |L4]], a new equation in quantum mechanics was introduced describing time evolution of this 

marginal distribution. The equation has completely classical form but contains all information 
about quantum system. 

The aim of this work is to consider two important types of nonclassical states of the trapped 
ion. First we consider squeezed and correlated states [|16| of the ion in the Paul trap || [F|. Then 
we study also even and odd coherent states of the ion irradiated by bichromatic laser field [p]]. 
The states of trapped ion are investigated in frame of symplectic tomography procedure suggested 
in |13|] and using the new quantum evolution equation of ||14j| . It is worthy to note that in |15 
endoscopy procedure to measure the quantum states of the trapped ion is elaborated in detail. 



1 0n leave from Lebedev Physical Institute, Moscow, Russia. 
Contribution to the Adriatico Research Conference "Interferometry 2" (ICTP, Trieste, 1-11 March, 1996), originally 
submitted in English 13 March, 1996. 



2 Gaussian Wigner Function 



The generic mixed squeezed state of the trapped ion with the density operator p is described 



by the Wigner function W (p, q) of the Gaussian form which contains five real parameters []17 



(the one-mode oscillator case). Two parameters are mean values of momentum (p) and position 
(q) and other three parameters are matrix elements of the real dispersion matrix m : 

win = 0>p ; 

m 12 = cr pq ; (1) 
m 2 2 = cr qq . 

Below, we will use the invariant parameters 

T = Tr m = a pp + a qq (2) 

and 

d = det m = a pp a qq - a 2 pq . (3) 

For one-mode system, the generic Gaussian Wigner function has the form (see, for example, IT71 
0) 

w (p, ti = -q ex P {~2d i aqqip ~ ^ 2 + app ( q ~ (g))2 _ 2cTpq ( p ~~ ( p) ^ q ~ ^] } • (4) 

The parameters (p) and (q) are given by the formulas 

(p) = Tr pp ; (5) 
(q) = Tr pg , (6) 

where the operators p and q are the quadrature components of creation a) and annihilation 
a operators 

a — a) 



a + a/ 



The matrix elements of the real symmetric dispersion matrix m are defined as follows 



(8) 



a pp = Tr pp 2 - (p) 2 ; 

a qq = Tipq 2 -(q) 2 - (9) 

Due to the physical meaning of the dispersions, the parameters a pp and a qq must be nonnegative 
numbers, so the invariant parameter T (|2]) is a positive number. Also the determinant d (0) of 
the dispersion matrix must be positive. For a pure Gaussian state, the parameter d = 1/4 . 



3 Squeezed States of Trapped Ion 

Since the ion in the Paul trap is described by the model of parametric oscillator in this section 
we review its properties. For the parametric oscillator with arbitrary time dependence of the 
frequency and the Hamiltonian 

where we put h — m — uj (0) = 1 and used expressions for the position and momentum operators 
in the coordinate representation, there exists the time-dependent integral of motion found in || 

A = ^=[s(t)p-i(t)q] , (11) 

where 

e (t) + uj 2 (t) e{t) = ; 

e(0) = 1; (12) 
e(0) = z, 

satisfying the commutation relation 

[A,A^} = 1. (13) 

For the trapped ion, the time- dependence of frequency is taken to be periodical one || 

uj 2 (t) = 1 + k 2 sin 2 Qt . (14) 

It is easy to show that packet solutions to the Schrodinger equation may be introduced and 
interpreted as coherent states 0, since they are eigenstates of the operator A (PD, of the form 

— i < T . . f Id 2 a 2 e*(t) V2ax) . . 

^:M) = *o(M)exp j-^--— ^ + , (15) 

where 

9 ( x ,t) = ^[e{t)r^e^^^ (16) 

is analog of the ground state of the oscillator and a is a complex number. Variances of the 
position and momentum of the parametric oscillator in the state (|16| ) are 

IfCQ I 2 . 

qq 2 ' 

(17) 

l^)| 2 



2 

and the correlation coefficient r of the position and momentum has the value corresponding to 



minimization of the Schrodinger uncertainty relation [|19] 



1 1 



&qq &pp 



4 1 - r 2 ' 

(18) 



a pq 



y/Pqq a pp 



If a qq < 1/2 (cr p p < 1/2) we have squeezing in quadrature components. 

The analogs of orthogonal and complete system of number states which are excited states of 
the ion in the Paul trap are obtained by expansion of (|15D into power series in a . We have 

*-<*■'> -Hr^^-Wi^n)- (19) 



and these squeezed and correlated number states are eigenstates of the invariant A^A . 
Another normalized solution to the Schrodinger equation 

*W(*,t) = 2tf«*oOM) ex p{-^r-i^} cosh T(?r' (20) 

where 

AT« = 6XP (l a |2 /2) , (21) 
2 \/cosh I a I 2 



is the squeezed even coherent state [Q] (the squeezed Schrodinger cat state). The odd coherent 
state of the parametric oscillator 

t( \, s , rM T / x ( I « I 2 e*(t)a 2 ) , y/2ax , . 

*Sr)(ar,t) = 2^-)*o(ar,*) expj-^J sinh , ( 22 ) 

where 

ATH = 6XP (l a |2 /2) , (23) 
2 Wsinh I a I 2 



satisfies the Schrodinger equation and is the eigenstate of the integral of motion A 2 (as well as 
the even coherent state) with the eigenvalue a 2 . These states are one-mode examples of squeezed 



and correlated Schrodinger cat states constructed in [20 



4 Tomography of Trapped Ion 

In [j!3| , it was shown that for the generic linear combination of quadratures which is a mea- 
surable observable ( h — 1 ) 

X = [xq + vp + 5, (24) 

where q and ft are the position and momentum, respectively, the marginal distribution 
w (X, fi, v, 5) (normalized with respect to the X variable), depending upon three extra real 
parameters //; v\ 5, is related to the state of the quantum system expressed in terms of its Wigner 
function W(q, p) as follows 

w (X, /i, v,S) = J exp [-ik(X - m -v V - 6)] W (q, p) . (25) 

As it follows from this formula, the marginal distribution depends on the difference of the variables 



So, we could introduce 



P (Y, //, v) = w (X, /i, v, 5) = w (X = Y, /i, u, 0) 



(26) 



which is marginal distribution depending on three variables. The physical meaning of the param- 
eters //; v\ 5 is that they describe ensemble of shifted, rotated and scaled reference frames in 
which the position X is measured. Formula (^5|) can be inverted and the Wigner function of the 
state can be expressed in terms of the marginal distribution |13| 

W(q, p) = (2tt) V exp (isX) w F {X, sq, sp, s) , (27) 

where wp (X, a, b, s) is the Fourier component of the marginal distribution taken with 
respect to the parameters //; v\ 5, namely, 



Wp (X, a, b, s) 



(2 



— - / w (X, /i, v, S) exp [— i(/ia + ub + 5s)] dp, dp d5 . 
7r) ,i J 



In [|14[], it was shown that for systems with the Hamiltonians 

i2 



H 



P 



+ V(q) 



(28) 



(29) 



the marginal distribution satisfies the quantum time-evolution equation 



d 

w — ii— — w — i 
ov 



V 



Id v d \ T , 

+ 1 - V 



d/ddd^i 2 88 



1 d v d 
d/d5dfl~ l 286 



w = . 



(30) 



The dot implies partial time derivative. This equation can be rewritten for the function P 
in the form 



V 



-1 d 



d/dY dfx 2 dY 



-Id v d 

+ i 77- 



\d/dYdfi 2dY / 



P = . 



(31) 



The measurable position Y is cyclic variable for the evolution equation. For the trapped ion, 
Eq. ( j30|) takes the form 

d 2 d 

w — fx — w + u (t)u — — w = , (32) 

OV OjJL 



or in terms of the function P, 



. 



(33) 



For squeezed and correlated states of the trapped ion, the Wigner function has the form (Q) in 
which quadrature variances and covariance for the state ( jl6|) are given by ( PD and ( ji8|) and 
quadrature means are equal to zero. For the state fll5D, the dispersion matrix is the same but two 
quadrature means have nonzero values 



<P> 



V2 

i 

71 



[as* + a*e) 



{as* + a*e) 



(34) 
(35) 



Calculating integral fl25|) one can show that for generic Gaussian packets of the trapped ion [also, 
for particular cases (|1^) and ([16]) ] the marginal distribution is 



1 f (X — X) 2 1 
w (X, /i, v, 5, t) = - exp \ — \ , (36) 

/2ir<T X (t) I 2 °x{t) J 



where the dispersion of the symplectic observable X and the mean value of the observable depend 
on time and parameters as follows 

a x (t) = fi 2 a qq + u 2 a pp + 2[iva pg ; (37) 
X = n(q) + u(p) + 5 . (38) 

It is the same form of Gaussian distribution discussed for free motion and harmonic oscillator 
in [[L4| but with different quadrature dispersions and means which are given by ([T?D , (|18"D and 
dH)) (HP ; correspondingly. 

One can check that the normalized marginal distribution (|36| ) with parameters ([J7|) and (|38]) 
and quadrature dispersions and means (Jl7|) , ([Tj|) and fl34]), (|35|) satisfies the evolution equation 



]). This follows from the observation that the marginal distribution depends on the symplectic 
parameters and time only through the dependence on the mean value X and dispersion ax{t) ■ 
Then calculating the first derivatives of the marginal distribution in all variables and using the 
Ehrenfest theorem for quadrature means and also the evolution equations for variances and co- 
variance of the parametric oscillator [JL^J , one can show that the evolution equation is satisfied by 
the function (|36|). 

The evolution of the Wigner function of trapped ion (as the system with quadratic Hamil- 
tonian) for any state is given by the following prescription (see, for example, []2T| ). Given the 
Wigner function W (p, q, t = 0) at the initial moment of time t = 0. Then the Wigner function 
at time t is obtained by the replacement 

W(p,q,t)= W[p(t),q(t),t = 0] , 

where time-dependent arguments are the linear integrals of motion of the quadratic system. The 
linear integrals of motion describe initial values of classical trajectories in the phase space of the 
system 



,^ e + e* e + e* 
P(t) = ~^—P ~ q] 



, . e — e e — e 
q{t) = — p + — q - 



(39) 



This ansatz follows from the statement that the density operator of the Hamiltonian system 
is the integral of motion and its matrix elements in any basis must depend on appropriate in- 
tegrals of motion. Consequently, for arbitrary mixed state of the ion in the Paul trap with the 
initial Gaussian Wigner function (£|), time evolution of the Wigner function is given by the same 
formula @j in which the variables p and q are replaced by the functions ([39D, correspondingly. 



The marginal distribution is related to the Wigner quasi-distribution function by (p5|). So, one 
can extend the above anzats to the marginal distribution. We formulate this anzats for the Paul 
trap as follows. Given an initial distribution function 

w (X, fi, v, 5, t = 0) = w (X - 5, fi, v) . (40) 

Then at time t ^ 0, the marginal distribution which is the solution to Eq. ( p2|) has the form 

w (X, /i, v,6,t)=wo(X- 5, fi(t), u(t) ) , (41) 

in which 



. . e — e e - e 
v (t) = v H u ; 



^ £ + £* £ + £* 
A 4 W = — n — ^ H ^ — M 



(42) 



5 Even and Odd Coherent States 

Now we will discuss marginal probability for nonclassical states of the parametric oscillator, 
such as even and odd coherent states 0. To describe such nonclassical states, as even and odd 
coherent states not only for an ion in the Paul trap but also for other types of trapped ions we 
consider construction of multimode even and odd coherent states (see, for example, |22| ). 



Symplectic tomography procedure [ 13 1 was extended for multimode systems in [E| 



We define multimode even and odd coherent states as 

| A±> = JVM ( | A)± | -A) ) , (43) 

where the multimode coherent state | A) is 

| A) =| en, or 2 , • • • , an) = D(A) | 0) , (44) 

and the multimode coherent state is created from the multimode vacuum state | 0) by the 
multimode displacement operator D(A). The definition of multimode even and odd coherent 
states is the obvious generalization of the one-mode even and odd coherent states. Normalization 
constants for multimode even and odd coherent states are 

N ( + ) = exp (| A | 2 /2) ; 

2 A /cosh I A I 2 



(45) 

expfl A | 2 /2) 
2 \/smh I A I 2 



where A = (ai, «2, • • • , ««) is a complex vector. 



The Wigner function for multimode coherent states is (see JT8| ) 



W A,B = 2N ex P ( ~ 2 ZZ * + 2 AZ * + 2 B * Z ~ AB * ~ ~ ~2) ' ^ 

where 

Z^. (47) 



For even and odd coherent states, the Wigner function is (see p2| ) 

W A± (q,p) = l^ (±) | 2 [^(A,B=A)( ( l'P) ± ^(A,B=-A)( ( l'P) 

± W (-A,B=A) (q» P) + W (-A,B=-A) (q» P) ] > ( 48 ) 
where explicit forms of are given by Eq. (f45[) . For multimode case, we use following notation 

AZ* = a x Z* + a 2 Z* 2 + ■■■ a n Z* ; 

(49) 

ZZ* = Z X Z{ + Z 2 Z* 2 + ■■■ + z n z* . 
If at time t = 0, one has the initial Wigner function of the system in the form 

W(p,q,* = 0) = W o (Q) , (50) 
the Wigner function of the system at time t is 

W(p,q,t) = W [A(t)Q + A(t)\. (51) 
where the matrix A (t) and the vector A (t) determine the linear integrals of motion for quadratic 



systems to which the parametric oscillator belongs |f2~T |. 

Thus, we get the expressions for the Wigner functions of even and odd coherent states of the 
trapped ion, if in formulas ( ^8|) we put n = 1 and replace p and q by the linear integrals of 
motion of the parametric oscillator 

A- At 

(52) 

A + At 
^2~' 



go 



where the invariant A is given by ([□]). The explicit form of the invariants po] qo which should 
be considered as c-numbers is given by (|39|) . Since the marginal distribution is related to density 
matrix which is the integral of motion, the time dependence of the marginal distribution for 
the parametric oscillator may be given by analogous prescription which was used for the Wigner 
function. If one knows the marginal distribution at the initial moment of time as the function of 



the parameters /i and u, at time t these parameters in the expression for the initial marginal 
distribution should be replaced by linear combinations of fi and v (|42|). 



The marginal distribution for initial even and odd coherent states was given in [|14 

1 



w ± {X, fi, u, 5) = iV (±)s 



7T 



(/X 2 + V 2 



{wi + w 2 ±w 3 ± w 4 } 



(53) 



where 



w 2 

w 3 = vol 

Y 

s 



cxp 



exp 



exp 



(Y 



/i 2 + v 2 
(Y + s + s* 



fi 2 + V 2 



-2 I a 



(Y — s + S* 
fi 2 + v 2 



\2l 



X-8 ; 

a/2 a (fi — iv) 



Thus, replacing the parameters v and fi in (|5^) by time-dependent functions ( f42"l) we get the 
marginal distribution of the trapped ion at any time moment t. 

The optical tomography procedure suggested in |I| is the partial case of symplectic to- 
mography, since for the optical tomography one has to take the symplectic parameters in the 
form 



;/ 
5 



COS ( 

sin^ 
. 



(54) 



Thus, the marginal distribution for rotated quadrature in case of initial even and odd coherent 
states takes the form 



w ± {X, (p, t = 
where 



N^ 2 — { Wl (X, <p)+w 2 (X, <p)±w 3 (X, <p)±w 4 (X, <p)} 



(55) 



Wi 




= exp 


- (X - 2 V2 | a 


COS (tp a - 


- V )f 








= exp 


- (X + 2 y/2 | a 


COS (ip a - 




> 






= exp 


-2 | a | 2 - (X - 


i 2 v2 sin 


(<Pa ~ 


< 




a 


= a | 


exp iip a . 









w 3 (X,p] 



This marginal distribution, used in the optical tomography scheme, at time t ^ is given by 
formula (|53|) in which the variables v and fi are replaced by the variables 



(56) 



e + e* e + e* . 
jj, (t, if) = — - — cos if H — sin if 



Using (|55| ) and ( |56"D and measuring the marginal distribution of homodyne output the Wigner 



function of the trapped ion may be reconstructed by means of the Radon transform [1C]. 

Thus, the rule formulated gives the time dependence of both the marginal distribution of 
rotated quadrature and the marginal distribution of squeezed and rotated quadrature for even 
and odd coherent states of trapped ion. Replacing /i and v in (KBI) by /i (t) and v (t) from 



(|42|), we get solution to the dynamical equation (|30"D for the ion in the Paul trap 



6 Conclusion 

In the present work, we calculated the marginal distribution of symplectic observable (which 
is generic linear quadrature) for nonclassical states of trapped ion modelled by the parametric 
quantum oscillator. Measurements of the marginal distribution give the possibility to measure 
the quantum states. In the case of particular choice of the parameters /i = cos<£>; v = sin<£>, the 
measurement is reduced to finding marginal distribution for homodyne output and reconstructing 



the Wigner function through the Radon transform of optical tomography scheme flU], [12| . The 
found distribution of generic linear quadrature satisfies the new classical-like equation of quantum 
dynamics introduced in the symplectic tomography scheme. 
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